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ABSTRACT 

We have investigated the properties of a bright K-band selected sample of early-type 
galaxies in the Hubble Deep Field, as representative of the field galaxy population. This 
dataset, based on public archives from HST and from deep observations at Hawaii and 
Kitt-Peak, is unique as for the morphological information on faint high-z sources, and 
for complete photometric and spectroscopic coverage. The majority of bright early-type 
galaxies in this field are found at redshifts z < 1.3 to share common properties with those 
of high-z cluster samples, as for the age and mass of the dominant stellar population - 
which are found to be as old as 3-5 Gyr and as massive as 10^^ Mq already at z ~ 1. The 
bright end of the E/SO population is already in place by that cosmic epoch, with space 
densities, masses and luminosities consistent with those of the local field population. 
There is no evidence of a steepening of the mass function from 2=0 to z=l, as inferred 
by some authors from analyses of optically-selected samples and favoured by hierarchical 
clustering models forming most of the E/SOs at 2; < 1. What distinguishes this sample is 
a remarkable absence of objects at 2 > 1.3, which would be expected as clearly detectable 
above the flux limits, given the aged properties of the lower redshift counterparts. So, 
something hide them at high redshifts. Merging could be an explanation, but it would 
require that already formed and old stellar systems would assemble on short timescales of 
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1 Gyr or less. Our analysis shows that, during such a major assembling, only a negligible 
fraction of the galactic mass in young stars has to be added, a rather contrived situation 
from the dynamical point of view. An alternative interpretation could be that a dust- 
polluted ISM obscures the first 3 to 4 Gyr of the evolution of these massive systems, 
after which a galactic wind and gas consumption makes them transparent. The presence 
of dust would have relevant (and testable) implications for the evaluation of the global 
energetics from galaxy formation and for the visibility of the early phases by current and 
future infrared facihties. 

Subject headings: galaxies: ellipticals - galaxies: ISM - dust, extinction infrared: galax- 
ies 
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1. INTRODUCTION 

Early-type galaxies have been studied in great de- 
tail in rieh clusters up to redsliift one and above (see 
Dickinson, 1997), with the basic result that passive 
evolution in luminosity and no evolution in mass is 
ruling them at least up to z = 1.2. Old ages and an 
early epoch of formation are then implied by these ob- 
servations. It is however still unclear how much repre- 
sentative of the general population is the rich-cluster 
environment: according to a common wisdom, early- 
type galaxies in the field might be much younger and 
follow a quite different evolutionary pattern (Baugh 
et al. 1996; Kaufmann et al. 1996). 

An exceedingly deep and clean view of field galaxy 
populations to high redshifts is provided by a long in- 
tegration of HST in the so-called Hubble Deep Field 
(HDF, Williams et al., 1996). Though having pro- 
vided important constraints on galaxy formation and 
evolution (Madau et al. 1997), the interpretation of 
these data is made difficult by the relatively short se- 
lection wavelengths (essentially U, B, V, I), which im- 
ply strong evolutionary and K-corrections as a func- 
tion of redshift and by the sensitivity to even small 
amounts of dust in the line-of-sight. The former two 
may be particularly severe for the early-type galaxies, 
because of their quickly evolving optical spectra and 
extreme K-corrections (e.g. Maoz, 1997). 

Major arguments for performing time-expensive 
surveys in the near-infrared, in particular in the K 
band, were to overcome most of these problems of 
optical selection. A further one is that at these wave- 
lengths the fluxes from all (low-mass) stars dominat- 
ing the baryonic content of a galaxy are counted. 

Various deep integrations have been performed in 
the HDF at near-IR (Cowie et al. 1996; Dickinson, 
1997) and mid-IR wavelengths (Rowan-Robinson et 
al. 1997; Aussel et al. 1997), which, combined with 
the extreme quality of the morphological informa- 
tion and the very good spectroscopic coverage, make 
this area unique, in particular for the investigation of 
early-type galaxies outside rich clusters. 

This letter is a preliminary report of an analysis of 
a bright complete K-band selected sub-sample of ellip- 
tical and SO galaxies in the HDF. Section 2 discusses 
physical properties of the galaxies as inferred from 
the morphology and from fits to the broad-band UV- 
optical-IR spectra. Section 3 investigates the statis- 
tical properties of the sample (counts, redshift distri- 
butions, identification statistics), and compare them 



with those of local early-type galaxies. The distribu- 
tion in space-time of the sample is found to display 
a discontinuity at ^; ~ 1.3, which apparently conflicts 
with the relaxed and passively behaving nature of the 
galaxies. Solutions of this paradox are discussed in 
Sections 4 and 5, by comparing the effects of merging 
with those of a long lasting dust-extinguished phase. 
We adopt Hq = 50 Km/ s/Mpc throughout the pa- 
per. 

2. THE FIELD POPULATION OF EARLY- 
TYPE GALAXIES 

2.1. The sample 

Our reference sample includes galaxies classified as 
ellipticals and SO in the Hawaii Active Catalogue of 
the HDF by Cowie et al. (1997), having K-band mag- 
nitudes brighter than 20 over the 4 square arcmin area 
covered by the 3 HST-WF cameras, observed in the J 
and H+K by Cowie et al. , and morphologically classi- 
fied by van den Berg et al. (1996). We have re-defined 
the HK-selected sample reported by Songaila et al. 
into a K-selected one by using the IRIM-KPNO im- 
age (Dickinson, 1997). Magnitudes are measured in 5 
arcsec diameter apertures. 

The complete sample consists of 19 early-type 
galaxies with K < 20, for 11 of which Cohen et al. 
(1996) report spectroscopic redshifts, while for the re- 
maining 8 a reliable estimate of z is obtained from fits 
of the optical-IR broad-band spectra. We have veri- 
fied that the typical error to be associated to the latter 
is not larger than 0.1. Morphologically, these galaxies 
appear as quite regular early-type systems up to the 
redshift limit of the sample, in no way dissimilar from 
local counterparts. 

The complete E/SO sample to if = 20 is supple- 
mented by a fainter one including 9 early-types and 
7 unclassified objects down to A" = 21. However, this 
enlarged sample has far less complete morphological 
and spectroscopic information, and is then used only 
to test the models for consistency at fainter magni- 
tudes and somewhat larger redshifts. 

2.2. Fitting the broad-band spectra 

The UV-optical-near /IR Spectral Energy Distribu- 
tions (SED) for all galaxies in the sample have been 
fitted with synthetic spectra based on a model briefly 
summarized here. The code follows two paths, one de- 
scribing the chemical evolution of the galaxy's ISM, 
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Fig. 1. Spectral energy distributions of four high-z early- type galaxies in the K-band HDF sample. Data on the 
source, the baryonic mass and the star-formation rate are indicated. 



the other associating to any galactic time a generation 

of stars with the appropriate mctallicity, and adding 
up the contributions light of all previous stellar pop- 
ulations. 

The chemical path adopts a Salpeter IMF with a 
lower limit M; = 0.15 Mq, and a Schmidt-type law 
for the Star Formation Rate (SFR): = vMgiff, 
where v is the SF efficiency. A further parameter 
is the time-scale tinfaii for the infall of the primor- 
dial gas. The library of isochrones is based on the 
Padova models, spanning a range in metallicities from 
Z = 0.008 to 0.05 (i.e. from 0.3 to 2.5 solar). The 
isochrones are modified to account for the contribu- 
tion by dust-embedded Mira and OH stars during the 
AGB phase (Bressan et al. 1997), when circum-stellar 
dust reddens the optical emission and produces an IR 
bump at 5-20 fim. This IR emission from AGB stars 
is important for galactic ages of 0.1 to a few Gyrs. 

A number of evolutionary patterns for the time- 
dependent SFR ^'(t) have been tried to reproduce the 
galaxy SEDs (four examples of which appear in Figure 



1). Though solutions are not univocally determined, a 

few general trends arc clearly written in the data. The 
first is that the bulk of the baryonic mass in objects 
observed at the typical redshift 2; ~ 1 has to be older 
than at least three Gyr, under the assumption of an 
high average stellar metallicity {Z ~ 2Zq). Lower 
metallicities would imply bluer spectra, hence even 
older ages. 

The second evidence is that, overimposed on this 
old stellar component dominating the spectra from 
0.6 to 2.2 fim and including typically more than 90% 
of the mass, a younger population is always present, 
as suggested by the often blue U and B colours. The 
evolution pattern (model 1 hereafter) fitting the SEDs 
predicts that a large fraction of the mass in stars has 
been generated during an intense burst lasting 0.7 
Gyr, after which the input of energy by stellar winds 
and supernovae produced a sudden outflow of the ISM 
through a galactic wind, stopping the bulk of the SF. 
The evolution at later epochs is mostly due to passive 
aging of already formed stellar populations. The dot- 
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Fig. 2. — The observed redshift distribution for 
sources identified as E/SO in the K-HDF sample (thick 
Une) versus predictions of model 1 (dotted line) and 
model 2 (shaded histogram). Most of the E/SO galax- 
ies are found at z > 0.5, because of K-correction and 
moderate evolution. 

ted lines in Fig. 1, fitting the red side of spectra, are 
the contribution of such an old population (though 
not exactly corresponding to model 1, but to a sim- 
ilar one described in Sect. 4 below). The global fit 
requires the additional contribution of a minor frac- 
tion of the mass in younger stars (dashed lines) . 

3. STATISTICAL ANALYSIS OF A K-SELECTED 
SAMPLE: THE PROBLEM OF THE RED- 
SHIFT DISTRIBUTION 

We have shown that a dominant fraction of the 

galactic baryonic mass has to be at least as old as 3 
Gyr for objects observed at 2 ~ 1. Correspondingly, 
the luminosity of these stellar populations is bound to 
increase with redshift during most of the cosmic time 
preceeding the epoch when galaxies are observed (the 
precise amount of such an increase depends weakly on 
the detailed values of the model parameters). 

It is of interest to check if this expected trend 
for the luminosity evolution can be reconciled with 
the statistical properties of this and other K-selected 
galaxy samples. To this end, model 1 has been com- 
bined with the the K-band luminosity function (LLF) 
of galaxies to predict sample statistics, under the as- 



sumption of a roughly constant galaxy mass function 
from redshift to 3. Our adopted LLF is taken from 
Gardner et al. (1997), which significantly updates 
previous determinations. The separate contributions 
of early- (E/SO) and late- type (Sp/Ir) galaxies have 
been defined using the morphological information on 
the optical LLF by Franceschini et al. (1988) and 
type-dependent B-K colours. 

We agree with Gardner et al. that the shape of 
the LLF can be naturally reconciled with faint K- 
band counts only in an open universe, while a clo- 
sure one would require some "ad hoc" evolutionary 
prescriptions. Then model 1 for E/SO galaxies (see 
Sect. 2.2), supplemented with a moderate evolution 
for Sp/Ir galaxies as in Mazzei et al. (1992), provides 
a good fit to the counts for go = 0.25. 

A crucial information is provided by the observed 
2;-distribution D(z) of early- type galaxies in the K- 
HDF (Figure 2). From ^ = to ^ = 1.2, D(z) is 
consistent with a moderate (essentially passive) evo- 
lution in an open universe (with 30=0.1-0.2). Again, a 
closure world model would require positive evolution, 
i.e. a larger comoving luminosity/density of massive 
E/SOs at z=l than locally. 

Hence, there is no evidence, up to this epoch, of 
an evolution of the baryon mass function. This is dif- 
ficult to reconcile with a decrease of the mass and 
luminosity due to progressive disappearence of big 
galaxies in favour of smaller mass units, as implied 
by some specifications of the hierarchical clustering 
scheme (e.g. Baugh, Cole, & Frenk, 1996). This also 
seem at odd with the finding by Kauffmann, Chariot 
& White (1996) of a strong negative evolution of the 
early-type population (by a factor 2-3 less in number 
density at z=l). 

Above z=1.2, however, model 1 has difficulties to 
reproduce the redshift distribution reported in Fig. 2 
and the identification statistics in Figure 3. A remark- 
able feature of the observational D(z) is the sharp cut- 
off at z > 1.3, such that no galaxies of the complete 
sample occur at such z. Model 1 is in no way able 
to explain the large number of sources at z = 1 fol- 
lowed by the rapid convergence above: the prediction 
would be of a much more gentle distribution, with less 
a pronounced peak and a large tail of galaxies (incud- 
ing typically half of sample objects) observable above 
z = 1.2. Lower values of qo would allow better fits 
of the observed D(z) up to ^; = 1.3, but would also 
worsen the mis- fit at higher z. 
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A similar effect is observed in the counts as a func- 
tion of the morphological class (Fig. 3). Here again 
sources identified as E/SOs show steep counts to K=19 
and a sudden convergence thereafter, while the model 
would predict much less rapid change in the slope. 
These problems remain forcing in various ways the 
model parameters. 




Fig. 3. — Identification statistics for E/SO (open 
squares) and Sp/Ir (open triangles) in the K-HDF, 
as a function of magnitude limit. Filled squares are 
the total observed counts. The dot-dashed line is the 
prediction for Sp/Ir, the dotted line is for E/SOs after 
model 1, while the long-dashed line for E/SOs after 
model 2. The thick continuous line is the prediction 
of model 2 for the total counts. 



4. EXPLAINING THE ENIGMA: DUST EX- 
TINTION DURING A PROLONGED AC- 
TIVITY PHASE 

4.1. Effects of merging at ^ > 1.3 

At redshifts larger than 1.3, early- type galaxies 
suddenly disappear from the K-selected sample, while 
we would expect to observe them to much higher 2, 
given the predicted luminosity enhancement when ap- 
proaching the most intense phase of star formation. 
The latter is expected to happen at z=2 to 5 after 
our analysis of the stellar populations in z=\ objects. 

In principle, this lack of high-z objects might be 
taken as an indication in favour of a fast merging- 



driven evolution of the mass function of galaxies at 
these redshifts. To be effective, however, this process 
has to operate on short timescales (< 1 Gyr) and, in 
particular, to assemble already formed, old and metal 
rich stellar populations, as implied by the spectra of 
2; ~ 1 — 1.2 galaxies (see Fig. 1). Only a small fraction 
(less than a few percent) of the mass in young stars 
has to be added to the galaxy, so stellar systems es- 
sentially deprived of gas are to be accreted, otherwise 
the gas would make stars actively during this dynam- 
ically perturbed phase. As suggested by N-body sim- 
ulations, typical outcomes of strong galaxy collisions 
are rings and bars, which are very rarely observed in 
HDF galaxies (van den Berg et al., 1996). 

Finally, it seems unlikely that the merge producing 
ellipticals and SOs should operate only among gas- 
poor systems, if we consider that gas-rich objects are 
so numerous in the field. 

4.2. Effects of dust during a prolonged SF 
pheise in field ellipticals 

An alternative interpretation, discussed here, for 
the lack of objects at z > l.'i in tlu^ K-HDF assumes 
that dust, present in a highly metal-enriched ISM dur- 
ing the first 3 to 4 Gyr of the evolution of massive 
galaxies, extinguishes the light emitted by high red- 
shift objects. The idea is that the SF efficiency in 
field galaxies with deep dark-matter potential wells 
is not so high to produce a galactic wind on short 
timescales after the onset of SF. Hence an ISM is 
present in the galaxy for an appreciable fraction of 
the Hubble time, during which the ISM is progres- 
sively enriched of metals and dust. While the gas 
fraction diminishes as stars are continuously formed 
at low rates, its metallicity increases and keeps the 
dust optical depht roughly constant with time. 

We have then extended the model of Sect. 2.2 to 
account for the effect of extinction and re-radiation 
by dust in star-forming molecular regions and more 
diffused in the galaxy body (see for details Granato 
et al. 1997). The amount of dust in the galaxy is 
assumed proportional to the residual gas fraction and 
chemical abundance of C and Si, while the relative 
amount of molecular to diffuse gas (Mm/Mg) is a 
model parameter. The molecular gas is divided into 
spherical clouds of assigned mass (5 10^ M©) and ra- 
dius (16 pc). Each generation of stars born within 
the cloud progressively escapes it on timescales of a 
few Myr with a velocity Vout ■ The emerging spectrum 
is obtained by solving the radiative transfer through 
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the cloud. Before escaping the gahixy, the hght aris- 
ing from young stars/molecular complexes, as well as 
from older stellar generations, further interacts with 
dust present in the diffuse gas component (the lat- 
ter dominates the global galactic extinction at late 
epochs, when the SF rate is low). The dimming of 
starlight and consequent dust emission are computed 
by describing the galaxy as a spherically symmetric 
system subdivided into volume elements, with radial 
dependences of stars and gas described by a King pro- 
file with a core radius Tc = 400 pc. 

This evolutionary code was used to explore al- 
ternatives to model 1. We found a solution (here- 
after model 2) with the following parameter values: 
baryonic mass M = IO^^Mq, (typical for our sam- 
ple galaxies), Unfaii = 0.1 Gyr, k = 1.5, v = 1.3, 
Mm/Mg = 0.3, Vout = 0.2 Km/ 8. Correspondingly, 
the SFR *(t) has a peak value of ~ 100 M^/yr at 
galactic time t = 0.2 Gyr and then decays as a power- 
law oc t{Gyr)~^ '^ . After 4.5 Gyr, the balance of gas 
thermal energy vs. input from supernavae breaks, and 
formally most of the small residual gas is then lost by 
the galaxy. After this event a low-level SF activity 
keeps on, due to either partial efficiency of the wind, 
or to stellar re-cycling. The average stellar metallicity 
of the remnant is Z = 2.5Zq. 

This scheme naturally accounts for the evidence 
(Fig. 1) that the highest redshift galaxies in the K- 
HDF show a combination of a massive amount of old 
(> 3 Gyr) stellar populations plus a small fraction 
of younger stars. Fits of galaxy SEDs from model 
2 (as reported in Fig. 1) can be found assuming an 
open world model with qo = 0.2. The redshifts zp for 
the onset of SF span a fairly wide range from 2 to 5 
(field early- types are not coeval in this sense, and in 
particular low- mass systems are younger on average), 
with an average value zp = 3.5. 

By constraction, the same scheme is successful in 
reproducing statistical distributions in Figs. 2 and 
3, in addition to other galaxy counts in the K-band. 
Now both the drastic cutoffs of D(z) for E/SOs at 
z=1.3 and of the counts at iiT > 19 are reproduced. 

Figure 4 summarizes our understanding of the 
E/SO sample galaxies, as for the baryonic mass (top), 
the residual SF rate ^ (mid), and the timescale 
needed by this SF activity to build the observed mass 
in stars (bottom). The residual SF rate is typically 
< 1 MQ/yr, and cannot contribute more than a few % 
of the baryonic mass for most of the galaxies. Note, 
finally, that the results in Fig. 4 are very robust: 




0.6 0.8 
redshift 

Fig. 4. — Basic properties of the E/SO sample as 
inferred from best-fitting model 2 on the observed 
SEDs. Top panel is an estimate of the total bary- 
onic mass; mid panel is the ongoing star formation 
rate; bottom panel is the ratio of the two, provid- 
ing an estimate of the time needed to synthesize that 
mass in stars with that rate of SF. 



essentially the same would be found using model 1 
instead of 2. 

5. CONCLUSIONS AND PERSPECTIVES 

Near-IR observations provide a view of galaxy evo- 
lution minimally biased by the effects of the evolu- 
tionary K-correction, dust extinction and changes in 
the M/L ratio due to the aging of stellar popula- 
tions. We have used a near-IR sample in the HDF 
with optimal morphological information and spectro- 
photometric coverage to study properties of distant 
early-type galaxies outside rich clusters. 

The SEDs of early-types in the HDF can be ex- 
plained in terms of a dominant old stellar population 
(older than 3-5 Gyr already at z=l) with a mass of 
10" Mq. The residual SF indicated by the blue U- 
B colours can only contribute some percent of the 
mass in stars. So, the main SF event producing the 
bulk of stars is confined to happen at ^; > 2 for rea- 
sonable world models. These general properties are 
not inconsistent with those found for cluster ellipti- 
cals studied at medium to high redshifts. Early-types 
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in the field are not coeval, however, the formation 
redshifts spanning a fairly wide range, zp — to 5, 
and the low-mass systems being younger on average. 

A basic problem is found when comparing the frac- 
tion of high-z objects expected on the basis of a simple 
photometric model with the observed distribution: a 
sizeable number would be expected at z > 1.2, while 
none is observed. Though not all redshifts in the sam- 
ple are spectroscopically confirmed, the latter conclu- 
sion seems quite robust: if optical spectroscopy has 
difficulties to enter this redshift domain, our anal- 
ysis, using broad-band galaxy SEDs when spectro- 
scopic redshifts are not available, has no obvious bi- 
ases, while the uncertainty in the photometric esti- 
mate of z for this kind of spectra is small. Of course, 
caution is in order about our conclusions because of 
the limited statistics. Improvements will be achieved 
by further areas deeply investigated by HSF (e.g. the 
southern HDF) and by further efforts of optical and 
near-IR spectroscopy. 
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Fig. 5.— Predicted spectra of a M = 3 10" M© 
galaxy at z = 2 observed at 0.5 (continuous), 1 (dot- 
ted), 4 (long dash), and 5 (short dash) Gyr, according 
to model 2. 

To explain the large number of E/SOs at z=l and 
their quick disappearence above, merging would have 
to operate on quite short timescales (< 1 Gyr) and 
to assemble mostly stellar systems, to avoid trigger- 
ing an excess amount of star formation. This is per- 



haps an uneasy situation as for the dynamics and for 
the probability of the event, that would be expected 
to favour, if any, dynamically efficient encounters be- 
tween the numerous gas-rich systems. 

We have discussed an alternative solution in terms 
of a somewhat protracted presence of a dusty and 
opaque ISM during the first 3-4 Gyr of galaxy's life, 
after which thermal unbalance in the ISM - due to a 
low residual gas mass and energy input by type-I SN - 
would clean up the object and make it transparent. A 
somewhat protracted SF activity in field, as opposed 
to cluster, ellipticals may be consistent with their ob- 
served "disky" morphology and the Hfj narrow-band 
indices. 

The implication of this is that long wavelength ob- 
servations would be better suited to detect and char- 
acterize the early evolutionary phases of spheroidal 
galaxies. Figure 5 illustrates possible spectra corre- 
sponding to different ages of a M = 3 10^^ Mq galaxy 
at z = 2, according to model 2. As indicated by the 
figure, various planned missions (in launch-time or- 
der, SIRTF, ESA's FIRST, NASA's NGST, and some 
ground-based observatories in exceptionally dry sites, 
like the South Pole) could discover unexpectedly in- 
tense activity of star-formation at z — 1.5 to 5. Also 
interesting tests of these ideas could be soon achieved 
with SCUBA on JCMT, and with ISO. 

This work has been funded by EC TMR Network 
Programmes FMRX-CT96-0068 and MRX-CT96-0086. 
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